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Abstract 



From a sample of 2722 ± 78 A+ decaying to the pK Tr"*" final state, we 
have observed, in the hadroproduction experiment E791 at Fermilab, 
143 ± 20 S° and 122 ± 18 S++ through their decays to A+7r=^. The mass 
difference M(S°)-M(A+) is measured to be (167.38±0.29±0.15) MeV; 
for M(S++) - M(A+), we find (167.76 ± 0.29 ± 0.15) MeV. The rate 
of A+ production from decays of the Sc triplet is (22 ± 2 ± 3) % of the 
total A^ production assuming equal rate of production from all three, as 
measured for S° and S^"*". We do not observe a statistically significant 
Sc baryon-antibaryon production asymmetry. The xp and pi spectra 
of A+ from Sc decays are observed to be similar to those for all A+'s 
produced. 
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Mass splittings among baryons containing different valence quarks provide a good 
probe of hadronic structure and of the forces that determine this structure. Isospin 
mass splittings have been attributed to the intrinsic difference between the light 
quark masses, and to the electromagnetic interactions among the quarks [1]. These 
splittings have been calculated with quark models and sum rules based on QCD [2-11]; 
in some of these calculations, it has been suggested that the color hyperfine mass 
splittings from the spin-spin interactions may also be important [12]. Experimentally, 
these mass splittings can be measured and different dynamical models can be tested. 
Especially interesting are the mass splittings in the charm baryon system, because the 
simplifications expected due to the large mass of the c quark allow more meaningful 
tests of models to be made. A disagreement of the Ec mass splitting measurements 
with sum rule predictions has recently been cited as posing a significant problem for 
the charm baryon sector of the quark model [13]. 

Measurement of charm baryon mass differences is most accessible through the 
hadronic decay of the 1=1 partner of the A^, Ec A^vr. Throughout this paper, the 
charge conjugate state is implicitly included. The Ec's are expected to be copiously 
produced in the hadroproduction of charm particles, and studies of the production 
characteristics of these charm baryon states provide information about the process 
of heavy quark hadronization. But to date, the only fixed-target hadroproduction 
experiment that has published results on Ec production is NA32 (ACCMOR) [14], 
based on a sample of 150 A^ pK~ir'^ candidates. They reported that from 230 GeV 
ir~ — N interactions, less than 5.3 % (90 % C.L.) of A^ decays come from E° and less 
than 5.2% (90% C.L.) from E^"*" decays. This can be compared to e"'"e~ [15-17] and 
photoproduction experiments [18,19] which report that the Ec triplet accounts for 
about 25 % of the total A^ production. 

In this paper, we report the observation of two of the three isospin states of the Ec, 
the E° and the E j""*" , and a measurement of the mass difference between each of these 

c c J 

two states and the A^. The results are based on 2722 it 78 A^ pK~ir'^ candidates 
produced in 500 GeV ir~ — N interactions in experiment E791 at Fermilab. Previous 
mass difference measurements have either been of low statistics or used several decay 
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modes of the added together [15-17,19], with different systematic uncertainties 
for each mode. We report measurements of the rate of production from Ec decays 
relative to the inclusive rate for A^. We also determine this relative rate as a function 
of Feynman-x [xp = Pz/p^""^) and transverse momentum squared (p^). 

The results given in this paper are based on the full dataset accumulated by 
experiment E791 in the 1991/92 Fermilab fixed-target run. The segmented target 
used in the experiment consisted of a 0.5 mm thick platinum foil followed by four 
diamond foils, each 1.6 mm thick. Each target center was separated from the next 
one by about 1.5 cm, thus allowing the observation of the decay of the A^ in air, 
minimizing the background from secondary interactions. The spectrometer [20] used 
in E791 was essentially the same as that used in experiments E691 and E769. It is 
a large-acceptance two-magnet spectrometer with eight planes of multiwire propor- 
tional chambers (MWPC) and six planes of silicon microstrip detectors (SMD) before 
the target for beam tracking, a 17-plane SMD system downstream of the target for 
track and vertex reconstruction, 35 drift chamber planes, two MWPC's, two multi- 
cell threshold Cherenkov counters, electromagnetic and hadronic calorimeters, and 
a muon detector. An important element of the experiment was its extremely fast 
data acquisition system [21] which was combined with a very open transverse-energy 
trigger to record a data sample of 20 billion interactions. 

The analysis begins with the selection of A^ pK~ir'^ candidate decays. For 
this, all events were first reconstructed and a candidate vertex position for the pri- 
mary interaction was determined. Events were selected for further analysis if at 
least one secondary vertex was also reconstructed. Tracks positively identified by the 
Cherenkov counters as kaons and protons were combined with the remaining tracks 
(assumed to be pions) in the event to search for A^ pK~ir'^ candidate decays. 

The final A^ sample was selected by the following requirements. The significance 
of separation of the candidate decay vertex from the primary vertex in the longitu- 
dinal direction was required to be greater than 7 o"a where a a is the measurement 
error on the longitudinal separation of the two vertices. The decay track candidates 
were each required to be inconsistent with having originated in the primary vertex 
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by requiring the transverse separation from the primary vertex to be greater than 
3o"f. The distance of closest approach between the hne of flight of the reconstructed 
decaying particle and the primary vertex was required to be less than 40 fim. The 
scalar sum of of the decay tracks, with pt measured relative to the direction of 
the parent, was required to be greater than 0.25 GeV^. We calculated, for each decay 
track, the ratio of its impact parameters to the secondary and primary vertices. The 
product of these ratios from the three decay tracks was required to be less than 0.01. 
Finally, the decay vertex had to be located at least 4 (where is the error on 
the position of the decay vertex along the beam direction) outside the target foils 
to remove secondary interactions. The flnal pK~ir'^ invariant mass distribution is 
shown in Fig. 1. From a flt to a Gaussian shape for the signal and a linear term 
for the background, a signal of 2722 it 78 events was obtained. Due to the low 
Cherenkov misidentiflcation rate and the good mass resolution of the spectrometer, 
the contamination from misidentifled D'^ and decays in the signal is less than 
5%. 

To search for the decay Ec A^vr, we add a pion of either charge to the pK~ir~^ 
state. This "bachelor" pion is required to be inconsistent with a kaon or proton 
according to information from the Cherenkov counters and to be not associated with 
another decay vertex reconstructed in the event. The mass differences between the 
Ec and A^ candidates, AM = M[pK~ Tr~^Tr^) — M[pK~Tr~^), are shown in Figs. 2 and 
3 for E° and E^"*", respectively, for pK~ir'^ combinations within 2.5 o"m of the fltted 
A^ mass, where aM is the measured width of the A^ signal (lOMeV). 

The mass difference spectra are flt to Gaussian distributions of flxed width and 
background shapes. The latter are relativistic phase space approximations, propor- 
tional to [1 + a{AM - m^)AM^] where a and /3 are free parameters in the flts. The 
width of the Gaussian signal term is determined by Monte Carlo studies [22] of the 
experimental resolution and flxed at 2.0 MeV for both the E° and E^"*". The result 
of an unbinned maximum likelihood flt to the mass difference plot for candidate E° 
decays (Fig. 2) yields 143±20 E° events with a mean AM of (167.44±0.29) MeV. For 
Fig. 3, the flt gives 122 ± 18 E++ events with a mass difference of (167.82 ± 0.29) MeV. 
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The natural widths of the E° and E^"*" resonances [12] are predicted to be small com- 
pared to the experimental resolution and are not included in the fit. Increasing or 
reducing the width of the Gaussian by 0.5 MeV changes the central value of AM by 
less than 0.1 MeV. 

We have studied the systematic uncertainty due to the background determination 
by three methods: 1) using different background parametrizations, 2) fixing the con- 
stants a and j3 to values determined by fitting the background events in the wings of 
the signal, and 3) using random combinatoric background given by a candidate 
from one event and a "bachelor" pion candidate from another event to determine a 
and j3. We found that the central value of the peak varies by less than 0.03 MeV. 
Since we measure the mass differences between the Ec and the A^ baryons instead of 
measuring the Ec mass directly, most of the experimental uncertainties from the re- 
construction of the A^ cancel. The remaining systematic uncertainty comes from the 
measurement of the momentum and decay angle of the "bachelor" pion. We use our 
measurement of the mass difference M(Z)*"'") — from D*'^ D^ir'^ decays to 

estimate this systematic error, since the momentum spectrum of the bachelor pions is 
similar and the mass difference technique is identical. Our measured mass difference 
M{D*+)-M{D°) is (145.47±0.01) MeV to be compared to the PDG fitted value [23] 
of (145.42 lb 0.05) MeV. Using our Monte Carlo, we find that the reconstructed mass 
difference for the D*+ - D° is shifted by (+0.05 ± 0.01) MeV relative to the generated 
mass difference. Similarly, comparing the reconstructed and generated values, the 
average shift in the mass difference between the Ec's and A^ is (+0.06 ± 0.04) MeV. 
We use (j) K'^ K~ decays (which have a similar Q to Ec decays) to estimate the 
overall uncertainty in our mass difference scale to be ±0.11 MeV. Finally, we estimate 
a systematic uncertainty of 0.07 MeV from the observed spread of the measured mass 
differences when the analysis cuts are varied. 

After correcting for the shift in the mass according to the Monte Carlo, our final 
result for the mass difference between E° and A+ is (167.38 ± 0.29 ± 0.15) MeV; for 
M(E++) - M(A+), the final result is (167.76 ± 0.29 ± 0.15) MeV. The corresponding 
values from the PDG fit are (167.3 ± 0.4) MeV and (168.04 ± 0.27) MeV. The only 
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recent results on these mass differences come from a photoproduction experiment 
which measured M(E°) - M(A+) to be (166.6 ± 0.5 ± 0.6) MeV and M(E++) - M(A+) 
to be (167.6ib0.6ib0.6) MeV [19]. Our results provide the best single measurements of 
these mass differences to date. We find AM (E++ -E°) to be (0.38 ±0.40 ±0.15) MeV 
which is in good agreement with the latest PDG value of (0.7 ±0.4) MeV. Theoretical 
predictions for this mass difference range from —18.0 to ±6.5 MeV and are listed in 
Table 1. The best agreement with data is obtained by Chan [4] who assumed in the 
framework of an SU(8) model, that hadron mass differences result from intrinsic quark 
mass differences and two-body spin-spin interactions together with the Coulomb and 
magnetic interactions. 

To calculate the rate of production from Ec decays relative to that for total 
inclusive A^ production, we use the Monte Carlo simulation to determine the re- 
construction efficiencies of A^, E° and E^"*". Our acceptance covers the kinematic 
region in the xp range above —0.1. Within this xp range, the ratio of the efficiency 
for reconstructing E°'"'""'" A^vr"'"'", A^ pK~ir'^ to the reconstruction efficiency 
for inclusive A^ pK~ir~^ decays is (66.5 ± 1.4) % and is the same for particle and 
antiparticle, and for both Ec charge states. Furthermore, this ratio of acceptances is 
found to be independent of xp within statistical errors so that our result is indepen- 
dent of the shape of the generated xp distribution of the A^. Assuming the branching 
ratio 5r(Ec ^ A+vr) = 1, we find that (7.9 ± 1.1 ± 1.0) % of A+ comes from E° and 
(6.7 ± 1.0 ± 1.0) % from E^"*" decays. The first error is statistical and the second one 
is systematic. The estimate of the systematic error comes from varying the width 
used in the fit for the mass difference between the Ec and the A^, using different 
background parametrizations, varying the analysis cuts and from the uncertainty in 
the tracking efficiency. 

If we assume equal production of states within the Ec isotriplet and average over 
the observed rate of production of E° and E^"*", then cr^^/a^+ = (7.3 ± 0.8 ± 1.0) % 
for each mode, or (22 ± 2 ± 3) % for all Ec production relative to A^ production. 
This is significantly higher than the previous upper limits reported for a 230 GeV 
ir~ beam [14]. Our 500 GeV result is consistent with the 25% seen in photoproduc- 
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tion [18,19] and e~'~e~ experiments [15,16] (Table 2). We have also checked for Ec 
baryon-antibaryon production asymmetry. In vr" — TV interactions, an enhanced pro- 
duction of leading D~ mesons has been observed compared to that for non-leading 
D'^ mesons in the forward region [24-26]. We measure a production ratio for E°/E^ 
of 1.34 ± 0.39 and a production ratio for E++/E" of 0.98 ± 0.30. We do not observe 
a statistically significant leading particle effect. 

To measure the relative production rate as a function oixp, we extract the number 
of decays in six bins of x^r in the range —0.1 to 0.5. We then impose a ±2.5 o"m cut 
on the pif "vr"'" mass and plot the difference between M{jpK~ tt^'k^) and M{jpK~ tt^) 
for the same range oixp- Each distribution is then fit to give the number of Ec decays. 
Similarly, we divide the signal into different bins between and lOGeV^ We nave 
used the Monte Carlo to check the ratios of acceptances between all and those 
which come from Ec decays as functions oi xp and p^. Within statistics, both ratios 
of the two acceptances are flat. After correcting for acceptances, we check the rate of 
production of A^ from Ec decays relative to the production for all A^ as a function 
oi xp and p^. The relative rate is found to be constant, showing no dependence on 
xp pI. This implies that the xp and pi distributions for A^ from Ec decays are 
similar to those from all A^, and subsequently that the Ec production spectrum is 
similar to the Ac spectrum (since the decay vr is from a low Q decay). 

In summary, we report the flrst observation of E° and E^"*" in vr" — AT interactions. 
We measure the difference between the E° and A^ masses to be (167.38 it 0.29 it 
0.15) MeV and the difference between E++ and A+ massses to be (167.76 ± 0.29 ± 
0.15) MeV. We observe an isospin mass splitting between the E^"*" and E° states of 
(0.38 ± 0.40 ± 0.15) MeV. We determine the rate of A+ production from the decay 
of the Ec triplet assuming isospin symmetry to be (22 it 2 it 3) % of the total A^ 
production in 500 GeV vr"— N interactions. 
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TABLES 



Reference 


AM(E++ - E°) in MeV 


Capstick [1] 


1.4 


Lane [3] 


-6.0 


Chan [4] 


0.4 


Hwang [5] 


3.0 


Sinha [6] 


1.5 ± 0.2 


Itoh [7] 


6.5 


Kalman [8] 


-2.7 


Lichtenberg [9] 


3.4 


Deshpande [10] 


-3.0 to -18.0 


Richard [11] 


-3.0 or 2.0 


Kwong [12] 


2.6 



TABLE L Theoretical calculations for AM(S++ - These predictions can be com- 
pared to our result of (0.38 ± 0.40 ± 0.15) MeV. 



Bar yon 


e+e- [16] 


145 GeV 7 [18] 


220 GeV 7 [19] 


230 GeV vr- [14] 


500 GeV vr" (our result) 


S° 


(9 ± 3) % 


(13±4±2)% 


(7.8 ±2.1)% 


< 5.3%(90%C.L.) 


(7.9 ±1.5)% 


S++ 


(9 ± 3) % 


(5 ± 3 ± 2) % 


(6.7 ± 1.9)% 


< 5.2%(90%C.L.) 


(6.7 ±1.4)% 



TABLE II. Rate of Sc production relative to total inclusive Ac production for different 
beams. 
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FIG. 1. Invariant mass distribution for candidate A+ pK~Tr'^ decays. The solid line 
corresponds to a fit to a Gaussian distribution plus a linear background. The dashed line 
shows the background alone. 
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FIG. 2. Mass difference spectrum, M(pif-7r+7r-) - M{pK-'K+), for S° pif-Tr+Tr" 
candidates with pK~Tr'^ combinations within 2.5 a of the fitted Ac mass. The solid line is 
a fit to a Gaussian distribution plus a background shape described in the text. 
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FIG. 3. Mass difference spectrum, M(pif"7r+7r+) - M{pK~'K+), for S++ pK~'K+'K+ 
candidates with pK~Tr~^ combinations within 2.5 ct of the fitted Ac mass. The solid line is 
a fit to a Gaussian distribution plus a background shape described in the text. 



15 



